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ABSTRACT 

, The theory of massive star formation currently suffers from a scarce observational 

base of massive young stellar objects to compare with. In this paper, we present inter- 
ferometric ^^CO (1-0), "CO (1-0), C^O (1-0), and 2.6 mm continuum images of the 
infrared source IRAS 19520+2759 together with complementary single-dish observa- 
tions of CS (1-0), obtained with the 34m antenna DSS-54 at the Madrid Deep Space 
I Communications Complex, as well as archive images at different wavelengths. As a 

result from our work, IRAS 19520-1-2759, with a controversial nature in the past, is 
firmly established as a massive young stellar object associated with a strong and com- 
pact millimetre source and driving a collimated outflow. In addition, a second fainter 
millimetre source is discovered about 4 arcsec to the south, which is also driving an 
outflow. Furthermore, the two millimetre sources are associated with C^^O clumps 
, elongated perpendicularly to the outflows, which may be related to rotating toroids. 

■ The masses of gas and dust of the millimetre sources are estimated to be around 100 

and 50 Mq. MMl, the dominant source at all wavelengths, with a total luminosity of 
(1-2) X 10^ Lq at 9 kpc, is however not associated with 6 cm emission down to a rms 
noise level of 0.1 mJy. We propose that IRAS 19520-1-2759 could be an example of the 
recent theoretical prediction of 'bloated' or 'swollen' star, i. e., a massive young stellar 
object whose radius has increased due to effects of accretion at a high-mass accretion 
rate. 

; Key words: Stars: formation - ISM: individual objects: IRAS 19520-^2759 - ISM: 
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1 INTRODUCTION 

Although high-mass (> 8 Mq) stars are clue pieces in the 
Universe and galaxy evolution, many questions remain open 
concerning their formation process. The root of the prob- 
lem is that the Kelvin-Helmholtz timescale for a high-mass 
star is much shorter than the free-fall timescale of the na- 
tal core, and thus the star reaches the main-sequence while 
still accreting matter. The important amounts of UV ra- 
diation emitted by the nascent massive star should halt 
at some point further accretion onto the star, limiting its 
final mass to a certain value which is currently not well 
established. From a theoretical point of view, there are a 
number of studies showing that it is possible to form a 40- 
50 Mq star through disc-mediated accretion (e.g., Yorke & 
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Sonnhalter 2002; Krumholz et al. 2009), and even it seems 
possible to form stars of up to 140 Mq (e. g., Kuiper et al. 
2010, 2011). However, from an observational point of view, 
a clear picture of the different evolutionary stages that a 
massive star undergoes until reaching the main-sequence is 
lacking, mainly due to the fast evolution of massive stars. 
Beuther & Shepherd (2005) and Keto (2007) suggest an ob- 
servational sequence for massive star formation, where the 
0-type stars progenitors are B-type objects embedded in 
large amounts of dust and gas and accreting at large mass in- 
fall rates (10~*-10~^ Mq yr~^). In this scenario, the B-type 
object drives a collimated molecular outflow, and while pro- 
gressively moving from early-B to late-0 and flnally early- 
O, the opening angle of the outflow progressively increases. 
Thus, the earliest stages in the formation of an 0-type star 
should be characterized by a central B-type object, driving 
a collimated outflow, and embedded in an infalling enve- 
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Figure 1. Left: Colorscale: Spitzer/MIPS at 24 /^m showing tiie surroundings of 119520 (colorscale units in MJy sr~^). Blaclc contours: 
20 cm emission from Landecker et al. (1990) and Kothes et al. (2006) tracing the Supernova Remnant (SNR). Contours are 7.5, 7.7, 7.9, 
and 8.1 K. Note that the 24 ^m emission is extended at the expanding front of the SNR and that 119520, located about 4 arcmin to the 
south, is the most prominent source in the field. Right: Zoom-in on the 119520 region, with the colorscale (in arbitrary units) tracing 
emission at 3.4 /im from the Wide-field Infrared Survey Explorer (WISE, Wright et al. 2010) and the white contours (of 0.12, 0.16, 0.20, 
0.24, 0.28, and 0.32 Jy beam"'^) corresponding to the 1.1 mm Bolocam emission (Rosolowsky et al. 2010; Aguirre et al. 2011). The two 
black stars indicate the position of the two millimetre sources detected in this work (Section 13.11 1. The linear scale, indicated by the 
vertical segment within the panels, is computed adopting a distance of 9 kpc to 119520 (Section 1). 



lope massive enough to provide the B-type object with the 
required mass to become an 0-type object. Typically, the 
massive young stellar objects (YSOs) in the earliest evolu- 
tionary stages also present emission from complex organic 
molecules (the so-called hot molecular cores, e.g., Cesaroni 
et al. 2010) and a lack of strong centimetre emission (e.g., 
Beuther et al. 2002; Molinari et al. 1998; Mottram et al. 
2011). A possible reason for the lack of an HII region is 
the quenching and/or trapping of the ionized gas by the in- 
falling dense envelope (e.g., Walmsley 1995; Molinari et al. 
1998; Keto 2002, 2003), or, as recently proposed by Hoare 
& Franco (2007), an increase of stellar radius as accretion 
proceeds onto the star due to trapped entropy (the so-called 
'swollen' or 'bloated' stars), which produces a decrease of the 
stellar effective temperature and of the UV radiation emit- 
ted by the nascent star (e.g., Yorke & Bodenheimer 2008; 
Hosokawa, Yorke, & Omukai 2010). 

Thus, the study of outflow and accretion signatures in 
massive YSOs is essential to progress in the understanding 
of massive star formation. However, there are only a few 
O-type YSO candidates (i.e., with bolometric luminosities 
> 50000 ^0) associated with disc-like structures or toroids 
of dense gas and possibly undergoing infall (e.g., Zhang, 
Ho, & Ohashi 1998; Pestalozzi et al. 2004; Solhns & Ho 
2005, Sollins et al. 2005a; Beltran et al. 2006; Beuther et 
al. 2007a,b; Beuther & Walsh 2008; Franco- Hernandez et 
al. 2009; Cesaroni et al. 2011; Furuya et al. 2011; Qiu et 
al. 2012; Jimenez-Serra et al. 2012; Wang et al. 2012), and 
it is not clear if these massive YSOs harbor true single O- 



type objects or a multiple system of early-B objects (e.g., 
Naranjo- Romero et al. 2012). This cannot be easily eluci- 
dated because most of the aforementioned massive YSOs 
undergoing infall are not detected at optical wavelengths 
and their spectral type cannot be assessed. On the other 
hand, YSOs with luminosities > 50000 Lq and detected in 
the optical/NIR have been found associated with accreting 
disks through near-infrared observations (e.g., Jiang et al. 
2008; Davies et al. 2010; de Wit et al. 2010; Kraus et al. 
2010; Wheelwright et al. 2012), but their number is also 
scarce. Thus it is extremely necessary to enlarge the sam- 
ple of O-type YSOs undergoing infall and actively accreting, 
and simultaneously detected in the optical/NIR to charac- 
terize not only the infalling envelope but also the central 
object. The results obtained from the joint analysis of mul- 
tiwavelength data presented in this paper indicate that one 
new actively accreting O-type (luminosity ~ 10^ -^0) '^''■n" 
didate is IRAS 19520+2759 (R.A. (J2000)=19''54™05';9 and 
Dec. (J2000)=+28°07'41"). 

The status of IRAS 19520-^2759 (hereafter, 119520) has 
been controversial: it was originally suggested to be an 
OH/IR star based on its infrared colors (Habing & Olnon 
1984). However, Arquilla & Kwok (1987) examined the In- 
fraRed Astronomical Satellite (IRAS) images towards this 
region and found the far-infrared emission to be extended. 
This, together with its spectral energy distribution (nearly 
flat from NIR to 100 jj,m), led these authors to conclude 
that this object was likely a young star embedded in an 
extended molecular cloud. In addition, Arquilla & Kwok 
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(1987) present a single-dish ^'^CO spectrum with broad wing 
emission (FWZI ~ 30 km s~^), similar to the hne wings as- 
sociated with outflows in star-forming regions, and Gledhill 
(2005) reports the detection of a polarized reflection neb- 
ula in the near- infrared. H2O and OH maser emission with 
single-peaked profiles (i.e., unlike double-peaked profiles in 
evolved stars) lends additional credence to the identifica- 
tion of this source as a YSO (Engels et al. 1984, Lewis et 
al. 1985). Moreover, the non detection of SiO maser emis- 
sion (Nakashima & Deguchi 2003) is also consistent with the 
YSO nature (e.g., Zapata et al. 2009b). 

119520 stands out in the infrared as the brightest and 
reddest object within a field of view of ~ 50 arcmin. A recent 
image from the Spitzer MIPSGAL survey (Multiband Imag- 
ing Photometer for Spitzer Inner Galactic Plane Survey, 
Carey et al. 2009) at 24 fim shows a very bright, compact 
source coincident with the IRAS source and surrounded by a 
much more extended (^15'), diffuse emission tracing a large 
molecular cloud (see Fig. [T|. Cold dust emission from this 
large molecular cloud is traced by Bolocam observations at 
1.1mm (Rosolowsky et al. 2010; Aguirre et al. 2011), which 
reveal an extended (~ 40 arcsec) source, of about 0.7 Jy of 
flux density, where 119520 is embedded (Fig. [l]). Interest- 
ingly, 119520 is located at the boundary of the supernova 
remnant G 65.1+0.6 in a region where the remnant is prob- 
ably interacting with the surrounding interstellar medium 
(Fig.[TJ Landecker et al. 1990; Tian & Leathy 2006). 

Hrivnak et al. (1985) identified an optical counter- 
part candidate of 119520, which appears as a bright, point- 
like source as observed with the Hubble Space Telescope 
(GO 9463, PI: R. Sahai). Recently, Sanchez Contreras et al. 
(2008) report echelle long-slit 3900-10900 A spectroscopy 
of 11952C1j, and show that the optical spectrum is domi- 
nated by emission lines consistent with the presence of a 
compact, photoionized region around an early-type (B or 
ealier) star (an in depth study of the optical spectrum is in 
preparation by Sanchez Contreras et al.). Finally, the intense 
Ha emission detected toward the nebular core shows broad 
wings (~2600 km s~^) and a remarkable P-Cygni profile 
with two blue-shifted absorption features, suggesting multi- 
ple stellar winds with (projected) expansion velocities of up 
to 800 km s~\ 

As already noted by Arquilla & Kwok (1987), the galac- 
tic longitude and negative Vlsr of 119520 indicate that this 
object is very remote. We estimate a value for its kine- 
matic distance of dk 8.8 kpc, determined from the tar- 
get radial velocity (Vlsr~— 16.5 km/s) and its galactic co- 
ordinates (Z=64.8131°, fe= -1-00. 1744°) by assuming a simple 
galactic rotation law and adopting a value for the A Oort 
constant of 14.4 km s~^kpc~^ and a galactocentric radius 
of 8.5 kpc (Kerr & Lynden-BeU, 1986). As indicated by Ar- 
quilla & Kwok (1987), this value of the distance is consistent 
with the non-kinematic distances, 8-10 kpc, of two Hll re- 
gions in the catalog of Blitz et al. (1982) that are likely 
to occupy approximately the same region of the Galaxy 
(given their similar I and radial velocities). The range of 
radial velocities displayed by the HI 21 cm-line emission 
profile, Vlsr=[— 20:— 26 km/s], at the boundary of the SNR 



Table 1. Parameters of line and 2.6 mm continuum OVRO inter- 
ferometric maps 





Array 


Beam size 


P.A. 


rms'^/rms' 


Maps 


Config. 


("x") 


n 


(mjy beam" 


12CO(1-0) 


CEH 


4.91 X 4.34 


85.5 


17/63 


i^CO(l-O) 


CEH 


5.36 X 4.57 


72.7 


12/25 


Ci«O(l-0) 


CEH 


5.56 X 4.54 


76.7 


10/12 


2.6 mm cont. 


CEH 


4.53 X 3.60 


83.6 


0.5 


2.6 mm cont. 


H 


2.29 X 1.45 


81.8 


0.7 



^ As measured in 1 MHz-wide channels with no signal 
^ As measured in the central 1 MHz-wide channel where the CO 
brightness distribution peaks 



G65. 1-1-0.6, next to which 119520 and the surrounding molec- 
ular cloud are located (Fig. 1), consistently yield a kinematic 
distance for the SNR of 9-10 kpc (Landecker et al. 1990; Tian 
& Leahy 2006). In view of these results, we adopt a distance 
of 9 kpc for 119520. 

The rest of this paper is organized as follows. In Sec- 
tion 2, we describe the observations, in Section 3 we give 
the main results from the OVRO 2.6 mm, and CO (1-0), 
"CO(l-O) and C^®O(l-0), in Section 4 we derive the main 
physical properties of the molecular gas, outflow and cen- 
tral object in 119520, and in Section 5 we discuss a possible 
interpretation for 119520 within the massive star formation 
framework. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 OVRO observations 

Interferometric mapping of the ^^CO, "CO, and C^O 
J=l— line emission at 108-115 GHz of 119520 was car- 
ried out using the six 10.4-m antennas millimetre array of 
the Owens Valley Radio Observatory (OVRO), which is now 
part of The Combined Array for Research in Millimeter- 
wave Astronomy (CARMA0). Observations were performed 
in different runs between 2002 and 2004 using the com- 
pact (C), equatorial (E), and high (H) array conflgura- 
tions, with baselines ranging between 4.01 and 93.3 kA in 
the uv plane. The coordinates of the tracking center are 
R.A. (J2000)=19''54™05':92 and Dec. (J2000)=-f28°07'41'.'5, 
and the total on-source time was '^18hr. 

Observations of the ^^CO, "CO, and C^®0 (J=l-0) 
transitions were performed simultaneously. The digital spec- 
tral line correlator was conflgured to provide a total band- 
width of 3x32 MHz ('^ 250 km s~^) with a channel spacing 
of IMHz (~ 2.6 km s"^) in the ^^CO and "CO lines. For 
C'^^O, the units of the cross-correlator were set to band- 
widths of 32 MHz (~85km s~^) with channel spacing of 
IMHz (2.7km s~^). The 2.6mm continuum emission was 
observed simultaneously using the dual-channel analog con- 
tinuum correlator. Although the continuum correlator pro- 
vides a total bandwidth of 4 GHz after combining both IF 
bands, our continuum maps have a bandwidth of 3 GHz since 



^ Due to its similar IRAS colors, 119520 was included in the list 
of pre Planetary Nebula candidates by these authors. 
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Table 2. Parameters of the sources detected with OVRO at 2.6 mm using the H configuration (highest angular resolution) 





Position^ 


Dec. ang. size 


Dec. P.A. 


fPcak a 


5, 


Mass^" 


Source 


a{J2000) 


(5(J2000) 


("x") 


(°) 


(mJy beam~^) 


(mJy) 


(Mq) 


MMl 


19:54:05.87 


28:07:40.9 


2.4 X 0.7 


-72 


9.0 ±0.7 


16 ± 5 


110 


MM2 


19:54:06.00 


28:07:36.3 


< 1.1 X 0.7 




4.8 ±0.7 


7± 2 


48 


MM2 W 


19:54:05.82 


28:07:36.3 


< 1.1 X 0.7 




3.4 ±0.7 


3± 1 


20 


All<= 


19:54:05.87 


28:07:41.1 


3.4 X 2.8 


52 


13.0 ±0.5 


20 ±6 


140 



^ Position, peak intensity and flux density are derived by fitting a Gaussian in the image domain (double gaussian for the case of MM2 
and MM2W). Uncertainty in the peak intensity is the rms noise of the cleaned image, a. Error in flux density has been calculated as 
yj' (ct Ssourcc/fbcam)^ + (""flux— scale )^ ) (Bcltran et al. 2001), where ^source and Sbcam S't'E the size of the source and the beam respectively, 
and cTflux-scalc is the error in the flux scale, which takes into account the uncertainty on the calibration applied to the flux density of 

the source {S„ X %uncertainty)- 

^ Masses derived assuming a dust temperature of 50 K (from a grey-body fit to the SED at wavelengths > 100 /im), and a dust (and gas) 
mass opacity coefficient of 0.003 cm^ g~^ (obtained by extrapolating the tabulated values of Ossenkopf &; Henning 1994, see Section [STTJ . 
The uncertainty in the masses due to the dust temperature and opacity law is estimated to be a factor of 4. 

Parameters derived from the image built using C+E+H configurations. For this case the peak intensity and fiux density are derived by 
computing statistics in a region of 12 arcsec of diameter. 



one of the 1 GHz bands, which contained the CO emission 
line, has not been used to generate the final maps. 

Data calibration was performed using the MMA soft- 
ware package (Scoville et al. 1993). Gain-calibration was 
done against a nearby quasar, J2025+337, which was ob- 
served at regular time intervals of 20 minutes before and 
after each target observation. The bright quasar 3C 454.3 
was observed at the beginning or end of the track for pass- 
band calibration. Absolute flux calibration was obtained by 
observing Uranus and 3C 454.3 as primary and secondary 
flux calibrators, respectively. Flux calibration errors could 
be of up to 20-30%, and absolute positional accuracies are 
estimated to be ~ 0.2 arcsec. 

Reconstruction of the maps from the visibilities was 
done using standard tasks of the Multichannel Image Recon- 
struction, Image Analysis and Display (MIRIAD) software 
(Sault et al. 1995). After Fourier transforming the measured 
visibilities with natural weighting, data were cleaned and 
maps restored. For the continuum, we have separately ob- 
tained two sets of maps, one set using all the data available 
from all three, C+E+H, configurations, and another one us- 
ing only data from the H configuration; the later provides 
maps with improved spatial resolution but lower S/N. For 
the line emission, imaging of the H configuration data suf- 
fers from an important decrease of S/N and missing flux, and 
for this reason we only show here the C+E+H configuration 
data for the line. The size and orientation of the synthe- 
sized clean beam together with the noise in our maps are 
given in Table [T] 



2.2 Single-dish observations with DSS-54 of 
Madrid DSCC 

We have used the DSS-54 antenna of NASA's Madrid Deep 
Space Communications Complex (DSCC) in Robledo de 
Chavela, under the 'Host Country' program. The antenna 
is a beam-wave-guide paraboloid, with a Coude focus, and 
has a diameter of 34 m. The observations were done using 
the dual-polarization Q-band High Electron Mobility Tran- 



sistor (HEMT) receiver attached to the antenna. We have 
observed both circular polarizations using the new wideband 
backend (Rizzo et al. 2012a). 

The CS(l-O) line (48.990955 GHz) was observed us- 
ing two spectral windows having bandwidths of 500 MHz 
and 16,384 channels, which provide a spectral resolution 
of 30.5kHz (0.19 kms"^ at the line frequency). Both 500- 
MHz bands were centred at 48.82 GHz in order to profit 
the wideband capabilities of the backend, and search for 
other molecular lines. Pointing was regularly checked, and 
was found to be better than 6 arcsec. The Half Power Beam 
Width (HPBW) of the antenna at the observed frequency 
is 44 arcsec. Observations were done in position-switching 
mode, with a reference located 3.6 arcmin in azimuth. 

Data were gathered on April 25, 2012, under good 
weather conditions; atmospheric opacity at the observed fre- 
quency was close to 0.1. The total on-source integration time 
was 46 minutes, and the system temperature was less than 
180 K. Estimated rms noise level was 40 mK, on a main- 
beam temperature scale. Control of the wideband backend 
and synchronization of the whole observing processes were 
done using the SDAI software (Rizzo et al. 2012b). The gath- 
ered data were written into FITS files by SDAI, and were 
subsequently read and further processed using the CLASS 
softwar^B 

2.3 VLA archive data 

We searched the vlaQ archive for centimetre continuum 
observations toward 119520. 119520 was observed at 6 cm 
(4.8 GHz) on 27 November 2003 with the array in the 

^ CLASS is part of the GILDAS software, developed by 
the Institute de Radioastronomie Millimetrique (IRAM). See 
rhttp://www.iram.fr/IRAMFR/GILDAS 

* The Very Large Array (VLA) is operated by the National Ra- 
dio Astronomy Observatory (NRAO), a facility of the National 
Science Foundation operated under cooperative agreement by As- 
sociated Universities, Inc. 
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Figure 2. Black contours: OVRO 2.6 mm continuum emission 
observed with the H configuration. Black contours are —3, 3, 
6, 9, and 12 times the rms of the map, 0.7 mjy beam~^. The 
synthesized beam of the H configuration data is shown in the 
bottom-right corner (2.29 X 1.45 arcsec^, P.A.= 81.8°). Grey con- 
tours: OVRO 2.6 mm continuum emission obtained by combining 
three configurations {C+E-(-H), providing a synthesized beam of 
4.53 X 3.60 arcsec^, P.A.=83.6°. Contours are -3, 3, 6, 12, and 
24 times the rms noise of the map, 0.5 mJy beam~^. The plus 
sign indicates the position of the Two Micron All Sky Survey 
(2MASS) point source. 
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Figure 3. Blue/red contours: OVRO ^^00(1-0) high veloc- 
ity emission. Blue contours correspond to integration in the 
velocity interval —50 to —26 km s~^, and red contours cor- 
respond to integration in the velocity —5.8 to 9.8 km s~^. 
Blue(red) contours start at 0.15 times the peak intensity, of 
13.02(8.18) Jy beam~^ km s~^, and increase in steps of 0.15 times 
the peak intensity. The synthesized beam of the ^■^CO (1-0) emis- 
sion, of 4.91 X 4.34 arcsec'^, P.A.=85.5°, is shown in the bottom- 
right corner. Green contours: VLA 6 cm emission. Contours are 
—3, 3, 6, and 9 times 0.11 mJy beam~^. The VLA synthesized 
beam is 3.84 X 3.54 arcsec^ with PA=88.8°. Black contours cor- 
respond to the 2.6 mm continuum emission as in Fig. [2] 



B configuration (project AH832, Urquhart et al. 200^, 
and at 2 cm (14.9 GHz) on 17 September 1986 with 
the array in the CnB configuration (project AP121). For 
project AH832, the phase center of the observations was 
RA(2000)=19''54'"05'!90 and Dec (J2000)=28°07'41'.'0. The 
absolute flux scale was set by observing the quasar 13314-305 
(3C286), for which we adopted a flux of 7.5 Jy at 6 cm. 
The quasar 2048-1-431, with a bootstrapped flux of 0.42 Jy, 
was observed regularly to calibrate the gains and phases. 
The total on-source time was ~2.5 minutes. For project 
AP121, the phase center was RA (2000)=19''54™05°40 and 
Dec (J2000)=28°07'37'.'6. The absolute flux scale was set by 
observing 3C286, with and adopted flux of 3.5 Jy at 2 cm; 
while the quasar 1923-1-210, with a bootstrapped flux of 
1.8 Jy, was used to calibrate the gains and phases. The 
source was observed for ~3.5 minutes. In both cases, the 
data reduction followed the VLA standard guidelines for cal- 
ibration of high-frequency data (or high-angular resolution 
data), using the NRAO package AIPS. Final images were 
produced with the robust parameter of Briggs (1995) set to 
5, corresponding to natural weighting, and applying a taper 
at 50 kA with the aim of recovering faint extended emission. 

^ The results of project AH832 are published in Urquhart et al. 
(2009), who report no detection in the field of 119520, with a 
rms noise of 0.16 mJy beam~^. However, in the cleaning process 
Urquhart et al. (2009) used a robust parameter of and did not 
apply any tapering to the data. Thus, with the aim of recovering 
possible faint and extended structure we re-did the calibration 
and imaging using a robust parameter of 5 (to improve the sen- 
sitivity) and tapering the data. 



The final la rms noise levels are 0.1 mJy beam~^ at 6 cm, 
and 0.6 mJy beam"'^ at 2 cm. The synthesized beams of the 
resulting images are 3'.'8 x 3'.'5 with PA=89° at 6 cm, and 
3'.'4 X 3'.'0 with PA=-82° at 2 cm. 



3 RESULTS 

3.1 OVRO 2.6 mm and VLA 6 cm continuum 

The 2.6 mm continuum emission toward 119520, after com- 
bining data from the 3 OVRO configurations, reveals one 
source in the field of view, which is slightly resolved, and 
whose peak position is coincident with the position of the 
infrared source (see Fig. [2] and Table [2}. Imaging using only 
the H configuration data (providing the highest angular res- 
olution) shows that the continuum emission splits up into 
at least two unresolved sources: one strong source detected 
at 14(7 (MMl) and one fainter source (MM2) about 4 arc- 
sec to the south detected at 6a, which has an elongation to 
the west (MM2W, Fig. [2]). The parameters of the millime- 
tre sources are listed in Table [2] For MMl, we measured a 
deconvolved size of 22000 x 6000 AU, at PA=-72°. As for 
MM2, a double-gaussian fit for MM2 and MM2W was re- 
quired to obtain a residual map with no excess of emission 
to the south of MMl, and both sources are unresolved. 

The total mass M of gas and dust from thermal contin- 
uum emission, assuming that the emission is optically thin, 
is: 
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Figure 4. OVRO ^^CO (1-0) channel maps for velocities ranging from —50 to 12 km s ^. Contours are —10, —5, 5, 10, 25, 50, 100, 150, 
200, 250, and 300 times the rms noise of the map, 0.017 Jy beam~^. The beam is shown in the bottom-right corner of the bottom-right 
panel, and the 2 stars mark the position of MMl and MM2. 



M 



B,{Ti)K,' 



(1) 



where 5*^ is the flux density at the frequency v, D is the 
distance to the Sun, Bv{Td) is the Planck function at the 
dust temperature Td, and is the absorption coefficient 
per unit of total (gas-|-dust) mass density. Writing Eq. (1) 
in practical units: 



M 



„0.048 i^/Td 



= 3.25 X 



- 1 







' d' 









(2) 



where Td is in K, u is in GHz, and is in cm^g~^. For the 
absorption coefficient at 115.27 GHz we used the extrapo- 
lated value from the tables of Ossenkopf & Henning (1994), 
for the case of thin ice mantles and density of 10® cm~'^, 
of 0.003 cm^g^^. Assuming Td of 50 K, the total mass 
estimated from the 2.6 mm continuum emission is around 



100 Mq for MMl, and 50, and 20 Mq for MM2 and MM2W, 
respectively. The global properties of the millimetre source 
(from the C-|-E-|-H configuration) are listed in the last row 
of Table [2] and, by adopting a mean molecular weight of 
2.8 mn, and a size (radius) of 13500 AU, we obatined an 
average density (of H2 molecules) of ~ 1.5 x 10® cm~"^, H2 
column density of ~ 4 x 10^^ cm~^, and A-^ ~ 440 ma^ 
(following Frerking, Langer, & Wilson 1982). 

Concerning the VLA centimetre emission, we only de- 
tected one unresolved source at 9cr about 40 arcsec to the 
southwest of MMl, as shown in Fig. [3l The flux density of 
this source (labeled 'VLAl') is ~ 1 mjy. At 2 cm we de- 



^ Note that among the 440 mag of visual extinction (derived 
from millimetre continuum data), only the half (~ 200 mag) are 
affecting the stellar light emitted towards the observer. 
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Figure 5. Top: Contours: zero-order moment (integrated inten- 
sity over all the velocity range) of ^^^00(1—0) towards 119520. 
Contour levels range from 10 to 100% of the peak intensity, 
37.07 Jy beam" ^ km s~^, increasing in steps of 10%. Bot- 
tom: Contours: zero-order moment of ^^CO(l-O) towards 119520. 
Contour levels range from 10 to 100% of the peak intensity, 
11.63 Jy beam~^ km s~^, increasing in steps of 10%. Colorscale: 
second-order moment (velocity dispersion, which has already been 
converted to linewidth) of ^^CO(l— 0). In both panels, star sym- 
bols correspond to the position of MMl and MM2 (Table (2)1 , and 
the synthesized beam is shown in the bottom-right corner. 



tected no source in the field shown in Fig. [S] preventing us 
from drawing any conclusion on the nature of VLAl. We 
also downloaded the most recent image at 6 cm of 119520, 
taken in 2008 within the framework of the CORNISbQ sur- 
vey (Hoare et al. 2012). In the 119520 region, the rms of 
the 6 cm image is 0.2 mjy beam~^, and we confirmed again 
the absence of any strong centimetre source associated with 
119520. 



3.2 OVRO ^^CO, ^=^CO and C^**0 emission 

Fig. [4] shows the channel map of the ^2 CO (1-0) emission, 
spanning 60 km s"'^ in total. The emission at high veloc- 
ities is compact, mainly found to the northeast or to the 
southwest of MMl, while the emission around the systemic 



http : / / www. ast .leeds .ac.uk /cornish / public / index. php 
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Figure 6. Top: White contours: zero-order moment (integrated 
intensity) of C^^O(l-O) towards 119520. Contour levels range 
from 20 to 100% of the peak intensity, 1.54 Jy beam~^ km s~^, 
increasing in steps of 10% (emission was integrated over all 
the velocity range). Colorscale: first-order moment (velocity) of 
C'^^O(l-O). Bottom: Idem as top with the colorscale showing 
the second-order moment (velocity dispersion, which has already 
been converted to linewidth). In both panels, black contours cor- 
respond to the H configuration 2.6 mm continuum emission as in 
Fig. [2] and the beam is shown in the bottom-right corner. 



velocity (—16.5 km s~^) is much more complex and extends 
about > 1 arcmin, being mainly elongated in the east-west 
direction. In particular, channels at —11.0 and —8.4 km s~^ 
show an X-shape structure to the west of MMl which could 
be tracing the cavity walls of an outflow. 

We computed the zero-order moment (integrated inten- 
sity) of the 1^00(1-0) emission for the entire velocity range 
(Fig. [5]-top), which reveals three main ^^CO clumps sur- 
rounded by extended emission: one associated with MMl, 
one about ~ 7 arcsec to the south, and the other > 10 arc- 
sec to east (there is also a possible fourth clump to the 
west of MMl). These three clumps are well separated in 
the zero-order moment of the ^^CO(l-O) emission, presented 
in Fig. [5]-bottom. The ^"^CO moment-zero map (integrated 
over all the channels where ^'^CO is detected, from —25 
to —11 km s~^) indicates that the strongest clump is the 
one associated with MMl, which was not obvious in the 
^■^CO map probably due to opacity effects. The other fainter 
^^CO clumps show however the broadest linewidths, of up 
to 10 km s"^ 
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Figure 7. CS (1-0) line observed with the DSS-54 antenna 
of the Madrid DSCC towards 119520. The spectrum has been 
smoothed to a spectral resolution of 0.38 km s~^. The dashed line 
indicates the velocity where the line peaks, —16.45 km s~^. 



We additionally computed the zero-order moment for 
^^CO in the high velocity range, integrating from —50 to 
—26 km for the blueshifted emission, and from —5.8 
to 9.8 km for the redshifted emission, and the result 
is presented in Fig. [S] The figure shows that the emission 
is bipolar and elongated roughly in the east-west direction, 
consistent with a bipolar flow driven by MMl. Note that to 
the south of MMl there is one blueshifted and one redshifted 
lobe to each side of MM2, suggesting that MM2 is driving 
also an outflow elongated in the southwest-northeast direc- 
tion. Since the clumps detected in ^''CO to the south, east 
and west of MMl span a broad velocity range, they could 
be associated with outflow emission as well. It is worth not- 
ing that the polarized reflection nebula detected by Gledhill 
(2005) is located about 2 arcsec to the east of 119520, where 
the blueshifted high- velocity ^^CO lobe is located, indicat- 
ing that possibly the CO bipolar outflow is creating a cavity 
which is seen in polarized light. This is found in other YSOs 
where outflows have excavated bipolar cavities (e.g.. Beck- 
ford et al. 2008). 

Finally, the C^**O(l-0) emission, detected in the same 
velocities approximately as ^^CO, is concentrated in one 
clump centred on MMl, with no clumpy emission in its sur- 
roundings. The zero-, flrst- (velocity field), and second-order 
(linewidth) moments for C^**0 are shown in Fig. |6l The in- 
tegrated C^®0 emission is elongated roughly in the north- 
south direction, with the southern part slightly tilted to the 
east, suggesting that we are observing the superposition of 
two cores (one north-south, and another one smaller and 
elongated in the southeast-northwest direction). It is inter- 
esting to note that MMl and MM2 lie at the center of these 
two C^^^O cores, and that the elongation of each of the two 
C^*0 cores is perpendicular to the corresponding outfiows. 
The velocity field of C^**0 reveals that MMl and MM2 are 
shifted by 2 km s^^. 



3.3 CS (1-0) with DSS-54 of Madrid DSCC 

The CS (1-0) line was detected using the DSS-54 antenna 
of Madrid DSCC (see Fig. [7)l , and the results of a Gaussian 
fit to the line are listed in Table [S] The linewidth of the line 
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Figure 8. Same as Fig. [T] (colorscale; 24 /^m emission; black 
contours: 20 cm emission) with the OVRO ^^CO high- velocity 
emission overlaid (red and blue contours). Note the bow shock 
structure 7 arcmin to the west of 119520, whose tail is pointing 
towards 119520. 



is ~ 2.3 km s~ , which is completely dominated by non- 
thermal motions, even if we assume that the dense core is at 
50 K. This broad linewidth could be produced either by tur- 
bulence injected by the outflow or by systemic motions such 
as rotation and/or infall. In fact, the line proflle is asymmet- 
ric, with the blueshifted side being slightly stronger than the 
redshifted side. This is the expected profile for an optically 
thick line tracing an infalling envelope (e.g., Torrelles et al. 
1995; Tsamis et al. 2008). 



3.4 Spitzer at 24 /im 

We downloaded the Spitzer (Infrared Array Camera (IRAC) 
+ MIPS) images of the 119520 field of view from the MIPS- 
GAL (Carey et al. 2009) databassi!]. In Fig. [S] we show a 
superposition of the '^^CO high- velocity OVRO emission on 
the 24 /im colorscale of Spitzer/MIPS. The MIPS image re- 
veals a structure with a bow-shock morphology about 7 ar- 
cmin to the west of 119520, whose tail is pointing towards 
119520, suggesting wind/jet activity in the east-west direc- 
tion as suggested also by the high- velocity '^^CO outflow. 
Such a wind/jet activity is consistent with the P-Cygni pro- 
file seen in the Ha line of 119520, from which a maximum 
expansion velocity between 400 and 800 km s~^ is inferred 
(Sanchez Contreras et al. 2008). Using an average velocity of 
600 km s~^, we estimate that the ejection event producing 
the bow shock took place about ~ 30000 tan(i) yr ago (with 
i the inclination with respect to the plane of the sky). 



http://irsa.ipac.caltech.edu/data/SPITZER/MIPSGAL/ 
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Figure 9. Top: Black (solid/dashed); l^co (1-0)/13CO(1-0) 
spectrum averaged over a region of 75 X 25 arcsec^ centred on 
MMl, used to derive tlie outflow parameters. Green: idem for 
MM2, for a region of 20 X 10 arcsec^. Blue and red vertical lines 
indicate where the outflow wing velocity range starts. Middle: 
Black (solid): ^^CO(l-O) averaged over a region of 20 arcsec (of 
diameter) centred on MMl. Red (long-dashed): idem for ^"^00(1- 
0). The vertical red dotted line indicates the velocity where the 
dip in 12 CO is seen. Bottom: Black (solid): 0^*0(1-0) averaged 
over a region of 15 X 10 arcsec^ centred on MMl. Green (dashed): 
01^0(1-0) averaged over a region of 10 X 7 arcsec'^ centred on 
MM2. Grey (solid): CS (1-0) spectrum as in Fig. [7| shown for 
comparison. In all panels, the vertical black dashed line indicates 
the systemic velocity. 



4 ANALYSIS 



4.1 



Physical parameters from ^^CO, ^=^00, C^'^O 
for MMl and MM2 



In Fig.[9]we show the ^^CO (1-0), ^=^00 (1-0), and C^^O (1- 
0) spectra averaged in different regions. The results from a 
Gaussian fit for all but the ^^CO transition are listed in 
Table [3] (corresponding to the spectra in the middle and 
bottom panels of Fig. [9]), together with the derived opacities, 
T, and column densities, A'^. Opacities for ^^CO and C^^O 
were estimated from equation: 



-In 



1 



where Tl is the intensity at the line center, Tex is the ex- 
citation temperature, Tbg is the background tamperature 
(adopted 2.7 K), and J^(r) = /ii//fe/(e'"'''*^ - 1) (with k 
being the Boltzmann constant, T the temperature and v the 
frequency). Tex, of around 9 K, was estimated from ^'^CO as- 
suming optically thick emission, and was found to be slightly 
smaller than the temperature obtained by Arquilla & Kwok 
(1987, 13 K). This difference is reasonable if we take into 
account the flux calibration uncertainties and the fact that 
Arquilla & Kwok (1987) observed with a single-dish antenna 
(with a HPBW of 45 arcsec) while we used the OVRO in- 
terferometer, being thus sensitive to different spatial scales 
(see below). 

We derived the ^^CO column density (from the 20 arc- 
sec averaged spectrum) following Palau et al. (2007), and 
using the opacity derived from "CO (1-0). The ^^CO (1-0) 
column density was estimated with the equation: 



isotp 



: 2.48 X 10" 



Tex 

IT 



Av 



km s 



1 _ ghi^/{fcTox) ' 



(4) 



(3) 



which can be used for 0^^0(1-0) as well. The obtained col- 
umn densities were around 2 x 10^^ cm"^ for ^^00(1-0), and 
(2-3) xlO^'' cm-2 for C^**O(l-0). As for the mass, we used 
the size listed in Table [3l adopted a mean molecular weight 
of 2.8, and abundances of XC^CO)^10-^, X("CO)=10-^ 
and X(Ci®0)=1.7 x 10"^ (Solomon, Sanders, & Scoville et 
al. 1979; Frerking, Langer, & Wilson 1982; Scoville et al. 
1986). The final masses are around 80 Mq including both 
MMl and MM2 (from "CO), and about 40 Mq and 20 Mq 
for MMl, and MM2, respectively (from C^^O, for which we 
could separate the emission from MMl and MM2). 

The systemic velocity is found at ~ —16.5 km s"^ 
(measured from ^^CO and CS), with MM2 being about 
2 km s~i redshifted with respect to MMl. The velocity is 
similar to that derived from H2O and OH maser (Engels et 
al. 1984; Lewis et al. 1985). The ^^CO spectrum (Fig. ^ 
middle) shows a double-peaked profile which could be due 
in part to the filtering of emission at systemic velocities 
and/or self-absorption of the cold foreground cloud. This 
is expected because ^^CO (1-0) is optically thick and traces 
gas of relatively low density (< 1000 cm~^). Ifowever, for 
opacity /filtering effects one would expect to see the maxi- 
mum absorption at the systemic velocity, while this is not 
observed. In fact, the ^^CO profile could be well explained 
if the material is infalling towards the central object, as the 
blueshifted ^^CO peak (and wing) is stronger than the red- 
shifted one, and the dip is redshifted by about 3 km s"^ with 
respect to the systemic velocity. Using the infall velocity de- 
rived from the ^^CO absorption dip, an approximate range 
of radii associated with this infall velocity (3000-36000 AU, 
from the measured size of the millimetre source and the 
C^^O clump. Tables [2] and [Sjl, and the H2 column density de- 
rived from i^CO (Table we estimated, following Beltran 
et al. (2006), a mass infall rate of (0.8-2.6) x 10"" Mq yr"^ 



4.1.1 Parameters of the CO outflows 

Table |4] lists the parameters estimated for the MMl and 
MM2 1^00(1-0) outfiows, computed following Palau et al. 
(2007; see also notes in the table; we used the ^^^CO spectra 
shown in Fig. [9]-top) and correcting for opacity effects. For 
the outflow inclination (with respect to the plane of the sky) , 
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Table 3. Physical parameters from averaged OVRO ^^CO (1-0), and C^^O (1-0) spectra, and DSS-54 single-pointing CS (1-0) 









Av 


Area 


A^isotp A^H2'' Sizc<: Mass 


Transition 




(K) (km s-i) 


(km s-i) 


(Kkms-i) t'^ 


(cm-2) (cm-2) (arcsec) (A/q) 


1^00(1-0) MM1+MM2 2.61 ±0.01 -16.54 ±0.02 


6.1 ± 0.1 


17.1 ±0.1 0.56 


1.8 X lO^"^ 1.8 X 


: 1022 11 75 


Ci8O(l-0) MMl 


0.56 ±0.02 -17.20 ±0.08 


5.5 ± 0.2 


3.3 ±0.1 0.10 


2.9 X 10^5 I J X 1022 8 36 


01*0(1-0) MM2 


0.38 ±0.03 -15.37 ±0.19 


5.5 ± 0.5 


2.2 ±0.2 0.06 


1.9 X 10^5 11 X 10^2 7 18 


OS (1-0) MM1+MM2 


0.26 ±0.03 -16.65 ±0.07 


2.3 ± 0.2 


0.64 ± 0.04 






Opacity estimated from the intensity at the line center 


, and using 


an excitation temperature of 9.3 K (derived from ^^CO for the 


optically thick 


case) . 












'° The adopted abundances are X(l''CO)=10 ® (Solomon, Sanders, & 


: Scoville 1979), and X(0l*0)=1.7 X 10" 


^ (Frerking, Langer, & 


Wilson 1982). 














FWHM. 














Table 4. Physical parameters of the outflows driven by MMl and MM2 








^dyn 


size N^2 Mout 


M ^ 


p a 


p> a 


pa 7" a 
-'^kin -^mcch 


Lobe 


(yr) 


(arcsec) (cm^^) (A^©) 


(Mq yr-i) 


(Mq km s-i) 1 


[Mq km s-i yr-l) 


(erg) (Lq) 


MMl-Red 


39000 


23 X 9 6.1 X 10^5 0.54 


1.4 X 10-5 


13 


3.3 X 10-* 


3.1 X 10*5 0.45 


MMl-Blue 


24000 


20 X 8 5.5 X 10^5 0.38 


1.5 X 10-5 


13 


5.5 X 10-* 


4.7 X 10*5 1.1 


MMl-AU 


32000 


1.2 X lO^** 0.92 


2.9 X 10-5 


26 


8.3 X 10-* 


7.8 X 10*5 1.6 


MM2-Red 


20000 


7x5 1.9 X lO^"^ 0.28 


1.4 X 10-5 


4 


1.8 X 10-* 


4.5 X 10** 0.11 


MM2-Blue 


9000 


4x4 3.0 X 10^'' 0.21 


2.3 X 10-5 


3 


2.9 X 10-* 


3.3 X 10** 0.14 


MM2-A11 


14000 


4.9 X lO^f* 0.49 


3.4 X 10-5 


6 


4.2 X 10-* 


7.8 X 10** 0.26 



^ Parameters are calculated foUowring Palau et al. (2007), for an inclination with respect to the plane of the sky equal to 45°, and are 
corrected for opacity using the ^^OO spectrum (see main text). For MMl we estimated an intensity in the line wings of ^^CO of about 
0.05 K (implying an opacity of ^^CO of 0.02 and a correction of ~ 2). From the line peak of the ^-^CO spectrum, ~ 3 K, we estimated 
and excitation temperature of ~ 6 K, assuming optically thick emission. For MM2 we estimated an intensity in the line wings of ^^OO 
of about 0.1 K (implying an opacity of ^^CO of 0.03 and a correction of ~ 3). From the line peak of the ^^CO spectrum, ~ 5 K, we 
estimated an excitation temperature of ~ 8 K, assuming optically thick emission. See Section l4.1.1l for further details. 



i, we adopted an intermediate angle of 45°, consistent with 
the bipolar outflow morphology. 

The dynamical times of the MMl and MM2 outflows 
are in the range 14000-32000 yr. In particular, the dynami- 
cal time of the outflow driven by MMl, of around 32000 yr, 
is very similar to the estimated timescale for the ejection 
event which produced the bow shock seen at 24 /im about 
7 arcmin to the west of 119520 (see Fig. [8]and Section [3.4(1 . 
reinforcing that 119520 probably had a recent and powerful 
ejection event. In fact, if the bow-shock is indeed driven by 
119520, this would be one of the longest jets known to date, 
extending for 18/cos(i) pc, similar to or even larger than 
the recently published updated length of the HH80-81 jet 
(Mcisque et al. 2012). As for the other outflow parameters, 
they are slightly smaller than the single-dish outflow param- 
eters derived by Arquilla & Kwok (1987, see their Table 2). 
This was expected because the interferometer filters out the 
most extended emission. For the case of our OVRO obser- 
vations, using the minimum uw-distance (see Section 2) and 
following Palau et al. (2010), we estimated that OVRO was 
filtering structures larger than ~ 23". From the CO (1-0) 
single-dish observations, Arquilla & Kwok (1987) estimate 
a total flux integrated over all the velocity range (single- 
point observations) of 1890 Jy km/s, and for the OVRO 
data we estimate (integrating both spectrally and spatially) 



870 Jy km s-\ Thus OVRO is recovering 46% of the single- 
dish flux, or missing a flux of 54%. 

In order to estimate the collimation factor of the outflow 
driven by MMl (i. e., excluding the outflow lobes of MM2), 
we followed the method used in Bally & Lada (1983) and 
Wu et al. (2004), where the collimation factor is the ratio 
of the major to minor radii of the outflow considered as 
an ellipsoid. By fitting a Gaussian to the redshifted lobe 
and another Gaussian to the blueshifted lobe, we obtained 
the minimum radius, of about 7 arcsec (deconvolved minor 
size of the gaussian fit), and the ma^dmum radius was es- 
timated from the total length of ~ 40 arcsec. This yields 
a collimation factor of 5.7, which is among the most high 
collimation factors measured in outflows driven by very lu- 
minous 10^ -^0) YSOs, either using single-dish telescopes 
(e.g., Wu et al. 2004; Qiu et al. 2011) or interferometer ar- 
rays (e.g., Cesaroni et al. 1999; Gibb et al. 2003; Beuther, 
Schilke, & Gueth 2004; Kumar, Tafalla & Bachiller 2004; 
Brooks et al. 2007; Zhang et al. 2007; Qiu et al. 2007; 2009, 
2012; Zapata et al. 2009a, 2011; Masque et al. 2012). 



4.2 Spectral Energy Distribution 

We built the Spectral Energy Distribution (SED) by compil- 
ing data from 2MASS (Skrutskie et al. 2006), WISE (Wright 
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Table 5. Results from the SED fit (Robitaille et al. 2007) 




X (|J,m) 

Figure 10. Observed and model SED of 119520 using Robitaille 
et al. (2007) tool. The filled circles show the observed fluxes. The 
thick coloured lines correspond to the best fit model (brown) and 
the three next best fits (grey). The models are calculated for an 
aperture of 48000 AU. The red thin line corresponds to the IRAS 
LRS spectrum. 



Parameter 


Best fit 
values 


Range of 
values^ 


Distance (kpc) 


8.7 


8.7-9.5 


A^'° (mag) 


4.4 


3.0-3.8 


M» (Mq) 


29 


26-33 


Tefl (K) 


39990 


38650-41500 




1.2 X 10^ 


(1.0-1.9) X 10^5 


Age (yr) 


61000 


(l.l-2.2)xl05 


Mdisc (Mq) 


0.3 


0.5-2 


Menv (Mq) 


138 


52-82 


Envelope outer radius (AU) 


10^ 


10^5 


Outflow opening angle dc (°) 


31 


13-27 


Outflow inclination (°)'^ 


59 


59-72 



^ Range of values for the three next best fits. 
^ Foreground interstellar extinction. 

Outfiow inclination with respect to the plane of the sky. 



et al. 2010), Akari (Murakami et al. 2007), IRAS (Neuge- 
bauer et al. 1984), Midcourse Space Experiment (MSX, 
Egan et al. 2003), Bolocam (Rosolowsky et al. 2010, Aguirre 
et al. 2011), and OVRO (this work). The global SED, shown 
in Fig. 1101 has a flat spectrum from 3 to 100 iim. In addition, 
we included optical photometry at 0.6 and 0.8 um extracted 
from Hubble Space Telescope images of I1952C|fj, and we also 
downloaded the IRAS Low Resolution Spectrometer (LRS) 
spectrc0 of 119520 (7.7-22.6 ^m). The LRS spectrum of 
119520 is classified by Kwok, Volk, & Bidelman (1997) as 
'unusual', indicating a flat continuum spectrum with un- 
usual features. The spectrum is quite noisy from 11 ^ra on, 
and shows no evidence of any silicate feature, apart from a 
very faint possible absorption at 9.7 ^.m, with a weak peak 
at shorter wavelengths. A weak silicate absorption with a 
weak peak at 8.5 ^m has also been found in other massive 
YSOs (e.g., Chen, Wang, & He 2000; CampbeU et al. 2006). 
New higher sensitivity observations should be carried out to 
properly study the spectral features of 119520 in this wave- 
length range. 

We have fitted the 119520 SED using the onUne tool 
provided by Robitaille et al. (2007) , which computes least- 
squares fits of pre-computed models of YSOs having discs 
and rotationally-fiattened infalling envelopes (with biconical 
outfiow cavities), to user-defined SEDs. We note that the 
parameters of the best-fit models of the online tool must be 

^ 119520 was imaged with the ACS/HRC instrument in the 
F606W and F814W broad band filters as part of a survey for 
candidate pre-planetary nebulae (GO 9463: Sahai et al. 2007), 
and appeared to be starlike in these images. 

IRAS Science Team (1986); for further deatils on LRS spectra, 
see http://irsa.ipac.caltcch.edu/IRASdocs/exp.sup/ to down- 
load the spectrum we used VizieR, Ochsenbein et al. (2000). 



regarded with caution, as recent work confronting the online 
tool results with spatial information data (e.g., Linz et al. 
2009; de Wit et al. 2010, 2011) indicate that the spatial 
information is required to break the degeneracy inherent to 
some of the online tool output parameters and to overcome 
the limited sampling of the parameter space (see Robitaille 
2008; Linz et al. 2009; Offner et al. 2012). Thus, we use the 
on-line SED fitting tool of Robitaille et al. (2007) as a first 
approach to the main properties of 119520. 

We set the input distance range to D = 8.5 — 9.5 kpc. 
The input interstellar extinction range was set to Av = 3 — 5, 
based on our estimate of A-^ ~ 4.1 ± 0.9 using the numerical 
algorithm provided by Hakkila et al. (1997), which computes 
the 3-dimensional visual interstellar extinction and its error 
from inputs of Galactic longitude, latitude, and distance, 
from a synthesis of several published studies. The best-fit 
model has D=8.7kpc, A^ — 4.4, a central star with mass 
M. = 29 Mq, effective temperature TcS = 39990 K, lumi- 
nosity L = 1.2 X 10® Lq, and age tagc ~ 6.1 x 10'' yr; a 
central disc mass Mdisc = 0.3 Mq, and an envelope mass 
Monv = 138 Mq. In this model, the (half) opening angle 
of the bipolar outflow cavity is 9c — 31°, and the outflow 
cavity axis is inclined at an angle of i = 59° to the plane of 
the sky. See Table [5] for the ranges of fltted values for the 
next three best-fit models. These models provide a good fit 
of the SED out to the far-infrared wavelengths. But at mil- 
limetre wavelengths, the models (based on a Monte Carlo 
code) do not have enough S/N, hence we are not able to 
make a comparison of the model and data for these wave- 
lengths. The circumstellar extinction is ^v(csm) = 2.4 in 
the best fit model, and lies in the range 3.1-4.1 for the next 
three best-fit models. This makes a total A^ (foreground -I- 
circumstellar) of ~ 7 mag for the best fit model. 
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5 DISCUSSION 



5.1 Comparison of the molecular gas and dust 

envelope properties of IRAS 19520+2759 with 
other massive YSOs 

In previous sections we have shown that the OVRO observa- 
tions reveal a strong dust condensation, referred to as MMl, 
with a total mass of around 100 Mq, and which is driving 
a high- velocity collimated ^^CO (1-0) outflow. The parame- 
ters of the outflow driven by MMl, if compared to those of 
outflows driven by high-mass YSOs, fit reasonably well the 
value expected for a massive YSO of ~ 10* Lq (e. g, Wu et 
al. 2004; Lopez-Sepulcre et al. 2009), after applying a factor 
of ~ 2 to account for the 54% of flux filtered out by OVRO 
(the works reporting correlations of outflow parameters use 
single-dish data and thus are sensitive to the most extended 
emission, which has been filtered out by OVRO). 

To further compare the 119520 envelope properties with 
previous works, we estimated the mass of the cloud associ- 
ated with 119520 by using the Bolocam image (Rosolowsky 
et al. 2010; Aguirre et al. 2011, see Fig. [T]), for which we 
measured a fiux density at 1.1 mm of 0.7 Jy, correspond- 
ing to 360-610 Mq, for a dust temperature in the range 
20-30 K and using the opacity law of Ossenkopf & Hen- 
ning (1994, ~ 0.0137 cm^ g"^ at 1.1 mm). This cloud mass 
and the CS (1-0) linewidth that we have measured (Table[3]) 
are consistent with the correlation found by Beuther et al. 
(2002) for a sample of high-mass protostellar objects. Sim- 
ilarly, the ^'^CO(l-O) linewidth that we have measured fol- 
lows the correlation with the bolometric luminosity found 
by Liu et al. (2010) for a sample of 98 massive YSOs. 

Finally, we plotted the values of the cloud mass and 
the measured bolometric luminosity on the mass-luminosity 
plot of Sridharan et al. (2002) and Molinari et al. (2008), 
and the values follow well the trend for massive YSOs. In 
particular, figure 9 of Molinari et al. (2008) indicates that 
119520 probably is in a transition stage between actively 
accreting objects and the 'envelope clean-up' phase. This 
is consistent with the fact that different methods used to 
estimate yield significantly different results, suggesting 
a non-homogenous density structure in the envelope and 
possibly the presence of holes and cavities. From the mil- 
limetre continuum emission we derived, assuming an homo- 
geneous density distribution, j4v~ 200 mag (Section 13. ip . 
This is about one order of magnitude larger than the 
estimated from the fit to the SED (total A^ is ~ 7, Sec- 
tion 14. 2p . and from the optical spectrum (ratio of Ha to 
H/? yields ~ 12 mag, Sanchez Contreras, in prep.), which 
are both along the line of sight. Overall, 119520 is proba- 
bly in a phase of main accretion and simultaneous start of 
envelope disruption as the jet interacts with the envelope 
creating a cavity. Examples of massive YSOs reported in 
the literature driving outfiows and with hints of simultane- 
ous envelope disruption are IRAS 05506-1-2414 (Sahai et al. 
2008) and V645 Cyg (Hamann & Persson 1989; Clarke et al. 
2006; Miroshnichenko et al. 2009), with 119520 being more 
luminous, and embedded in a one order of magnitude more 
massive envelope. 



5.2 The central star of IRAS 19520+2759 

The agreement between the position of the optical/infrared 
source and the millimetre source suggests that they are trac- 
ing the same object. The position of the 2MASS source is 
RA(J2000): 19:54:05.86, Dec (J2000): +28:07:40.6. The po- 
sition of MMl (Table[2ll is offset with respect to the 2MASS 
position by (+0.15, +0.27) arcsec, which falls well within 
the OVRO (~ 0.2 arcsec) and 2MASS (~ 0.1 arcsec, Skrut- 
skie et al. 2006) positional uncertainties. In addition, the 
preliminary analysis of the optical spectrum of 119520 pre- 
sented in Sanchez Contreras et al. (2008) points to a hot ob- 
ject photoionizing the circumstellar gas and driving a dense 
and compact wind. Thus, the properties of both the optical 
source and the millimetre source indicate that we are dealing 
with a young and massive object. 

Further hints on the nature of the optical source may 
come from the non-detection of centimetre emission: if the 
central star was a (main-sequence) 09 star with a typical 
lonizmg fiux iVi ~ 10** s"^ (e. g., Panagia 1973; Smith, Nor- 
ris & Crowther 2002), the expected 6 cm flux density from 
an optically thin HII region, with an electron temperature 
of ~ 10000 K and placed at a distance of 9 kpc, would be 
100 mjy, which should have been well detected with the 
VLA observations (Section 12. 3|F^ . If the spectral type is Bl 
(or later), the non-detection of the 6 cm emission is natu- 
rally explained, as the fiux of ionizing photons in this case 
is smaller than for the 09 case by one-two orders of magni- 
tude. 



5.3 IRAS 19520+2759: a B-type star becoming an 
O-type star? 

As discussed above, the centimetre properties of 119520 
are most consistent with an early-B type classification for 
the central star; however, the spectral type inferred from 
its high luminosity, ~ 10'' Lq, assuming it is a zero-age- 
main-sequence star, is significantly earlier, 06-07 (e.g., 
Bernasconi & Maeder 1996; Martins et al. 2005). The prop- 
erties of 119520 can be reconciled if the object is accreting, 
which is strongly suggested by several observational facts: i) 
the infrared source is associated with a compact and dusty 
envelope, MMl, of ~ 100 Mq; u) MMl is driving a high- 
velocity collimated outflow with a dynamical timescale of 
~ 30000 yr; and iii) the OVRO C^*O(l-0) emission, with a 
size of ^ 70000 AU, is elongated perpendicular to the out- 
flow, similar to the molecular toroids found in other mas- 
sive star-forming regions (e.g., Beltran et al. 2005, 2011; 
Beuther et al. 2007b; Zapata, Tang, & Leurini 2010; Furuya 
et al. 2011), which usually show a rotation velocity pattern. 



If we assume optically thick emission, the 6 cm upper limit can 
only be reproduced with an O-type star if the free-free emission 
comes from a region of very small size (<SC 1000 AU) and large elec- 
tron density (3> 10® cm~^), which is very unlikely because such 
a small size (typical of hyper-compact HII regions, e.g., Kurtz 
2005) would require a very dense envelope strongly infalling and 
quenching the ionized gas (e.g., Walmsley 1995; Molinari et al. 
1998), preventing the object to be detected in the optical (e.g., 
GlO.6-0.4: Sollins & Ho 2005; G28.20-0.05: SolUns et al. 2005; 
NGC7538-IRS1: Sandell et al. 2009). 
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All this suggests that 119520 has been recently undergoing 
active accretion. 

Hoare & Franco (2007) have proposed, and numerical 
simulations have shown (e.g., Yorke & Bodenheimer 2008; 
Hosokawa & Omukai 2009; Hosokawa, Yorke, & Omukai 
2010), that the main effect of accretion onto a massive star 
(with mass accretion rates up to 10~^ Mq yr~^), regardless 
of the geometry of the accretion, is an increase of the stel- 
lar radius, giving rise to the so-called 'swollen' or 'bloated' 
stars. The 'swollen-star' scenario naturally explains the low 
ionizing flux observed in 119520, even though with a bolo- 
metric luminosity of 10^ Lq . In this scenario, the stellar ra- 
dius attains its maximum value when the star has accreted 
~ 10 Mq (Hosokawa, Yorke, & Omukai 2010). Assuming 
that this is the current value of the stellar mass of 119520, 
and given the mass measured for the envelope of ~ 100 Mq 
(Section 3), the total mass which could still be accreted onto 
the star, for a star formation efficiency of about ~ 10-20%, 
is 10-20 Mq, implying a final stellar mass of 20-30 Mq, 
which corresponds to an O-type star. 

The 'swollen-star' scenario has also been proposed for 
other massive YSOs. Davies et al. (2011) require a 'swollen- 
star' phase to reproduce the observed luminosity distri- 
bution of massive YSOs in the sample of the Red MSX 
Source survey (RMS, e.g., Urquhart et al. 2008; see also 
Mottram et al. 2011); Simpson et al. (2012) and Bik et 
al. (2012) explain the low effective temperatures or offsets 
from main-sequence in the Hcrtzsprung Russell diagram 
by invoking a 'swollen-star' phase for the massive YSOs 
of G333.2— 0.4 or W3Main complexes. Finally, individual 
objects studied in detail within the 'swollen-star' scenario 
are B275 (Ochsendorf et al. 2011), M8E-IR (Linz et al. 
2009), CRL2136 (de Wit et al. 2011), and OrionKLlRc2 
(Morino et al. 1998; Testi, Tan, & Palla 2010). 119520, with 
a bolometric luminosity around 10® Lq , would be the most 
luminous of the four aforementioned 'swollen- star' candi- 
dates. Accurate spectral type classification of the central 
star through, e. g., optical spectroscopy in the 4000-5000 A 
range, is needed to firmly establish the 'swollen' star hy- 
pothesis proposed for 119520. 



6 CONCLUSIONS 

We present the results of OVRO observations at 2.6 mm and 
^^CO (1-0), "CO (1-0), and C^O (1-0) transitions, as well 
as observations of CS (1 0) with the 34 m antenna DSS-54 
of Madrid DSCC, towards IRAS 19520-1-2759, a very bright 
infrared object embedded in a massive cloud. Our main con- 
clusions can be summarized as follows: 

- The 2.6 mm continuum emission is dominated by one 
strong millimetre source, MMl, clearly associated with the 

bright infrared source, and which is only barely resolved, 
with a size of 22000 x 6000 AU, PA=-72°, and a mass of 
~ 100 Mq. About 4 arcsec to the south there is a fainter 
and unresolved source, MM2, of about ~ 50 Mq . 

- Both MMl and MM2 are driving high- velocity ^^CO 
outflows, with outflow parameters typical of massive YSOs, 
and the outflow driven by MMl shows a high collimation 
factor. In addition, both MMl and MM2 are associated with 
^^CO and C'^^O emission, whose structure, especially for 



C O, is elongated perpendicularly to the outflow directions, 
similar to toroids seen in other massive star-forming regions. 

- The CS line profile indicates that the object is em- 
bedded in dense gas with important contributions of non- 
thermal motions. 

- We fitted the SED of MMl using the SED fitting tool 
of Robitaille et al. (2007). The best fit yields a foreground 
extinction of 4.4 mag, an envelope mass of 140 Mq, and 
a bolometric luminosity of ~ 10® Lq. For this bolometric 
luminosity one would expect that the star should emit an 
ionizing flux high enough to detect the free-free emission at 
centimetre wavelengths. However, no source associated with 
IRAS 19520+2759 is detected at 6 cm, with a rms noise of 
0.1 mjy, which could be due to on-going accretion onto the 
star. 

Thus, we present evidence that IRAS 19520+2759 is 
a luminous YSO embedded in a massive and dense enve- 
lope, which is driving a coUimated outflow and does not 
show strong centimetre emission. We propose that this 
could be an example of a 'swollen' or 'bloated' star as pro- 
posed by Hosokawa, Yorke, & Omukai (2010), where accre- 
tion produces an increase of the stellar radius, and a de- 
crease of effective temperature and ionizing flux. Overall, 
IRAS 19520+2759 seems to be an excellent accreting mas- 
sive YSO candidate which could form an O-type star in the 
future. 
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